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Abstract 

We study bifurcations of a three-dimensional diffeomorphism, go, that has a quadratic 
homochnic tangency to a saddle- focus fixed point with multipliers (Ae*''', Ae"*"^, 7), where 
< A < 1 < I7I and |A^7| — 1. We show that in a three-parameter family, g^, of dif- 
feomorphisms close to go, there exist infinitely many open regions near e = where the 
corresponding normal form of the first return map to a neighborhood of a homoclinic point 
is a three-dimensional Henon-like map. This map possesses, in some parameter regions, a 
"wild-hyperbolic" Lorenz-type strange attractor. Thus, we show that this homoclinic bifur- 
cation leads to a strange attractor. We also discuss the place that these three-dimensional 
Henon maps occupy in the class of quadratic volume-preserving diffeomorphisms. 

1 Introduction 

In this paper we are concerned with the study of the three-dimensional map (x, y, z) = f{x, y, z) G 
defined by 

x = y, y = z, z = Mi + Bx + M2y - z^ , (1) 
with three parameters. Mi, M2, and B. This map is a natural generalization of the famous two- 



dimensional Henon map H76 ; indeed, like the latter, the map is quadratic, has constant 



Jacobian det(D/) = B and moreover, when B = 0, reduces to the two-dimensional Henon map. 
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schools" No. NSh-838. 2003.2. JDM would like to acknowledge support of the NSF through grant DMS-0202032. 
This collaboration was supported by the CRDF under grant No. RU-M1-2583-MO-04. 
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Therefore, it is natural to call the Henon map." We will study a number of properties 
of including: 

1. its nontrivial dynamics and, most importantly, the existence of wild-hyperbolic strange 
attractors; 

2. its origination in homoclinic bifurcations; 

3. its structure as a member of the family of three-dimensional quadratic maps with constant 
Jacobian. 

Three-dimensional maps have not been as widely studied as one- and two-dimensional maps, 
though in recent years there has been a considerable increase in their study. One of important 
reasons for this fact is that multi-dimensional dynamical systems (the dimension of the phase 
space is at least four for flows and three for maps) can exhibit complicated dynamics that is 
decidedly distinct from lower-dimensional cases. In particular they can possess a new variety of 
strange attractor called a wild-hyperbolic attractor |ST98j . 

Wild-hyperbolic attractors have no uniform hyperbolic structure; moreover, they may con- 
tain homoclinic tangencies and coexisting periodic orbits of different types. However, these 
attractors contain no stable invariant subsets such as attracting periodic orbits. Most crucially, 
the latter property holds for open regions in the space of smooth dynamical systems. Thus, 
wild-hyperbolic attractors are very interesting objects and their detection in concrete models is 
a fundamental problem of nonlinear dynamics. 

Most of the strange attractors occurring in two-dimensional, smooth maps (or three-di- 
mensional flows) are not structurally stable since they possess homoclinic tangencies. For such 
systems, even though the orbits are observed to be quite chaotic for some parameter values, very 
small changes in the parameters may destroy the aperiodicity and give rise to periodic behavior. 
Such parameter regions are called stability windows. Even if these windows are not "observable" 
in numerical computations, their appearance (under arbitrarily small perturbations) is inevitable 
because of the prevalence of homoclinic tangencies. Such attractors were called quasi- attractors 
in |AS83j . One of the most challenging problems in dynamical systems is to distinguish aperiodic 
motion from extremely long periodic motion. Moreover, since most of the chaotic attractors we 
see in applications are definitely quasi- attractors, it is unclear in which sense they are chaotic, 
or how one can define and measure the "probability" for a quasi-attractor to be chaotic. 

However, it is well-known that there are strange attractors that are free of these problems; for 
example, hyperbolic attractors and the Lorenz attractor. Though the latter is not structurally 
stable |GW79j . |ABS82j (in contrast to hyperbolic attractors) it contains no stable periodic orbits, 
every orbit in it has positive maximal Lyapunov exponent, etc. Moreover, these properties are 
robust |ABS77j . |GW79j . The reason is that Lorenz attractor possesses a pseudo- hyperbolic 
structure. In terms of Turaev and Shilnikov 'ST98^ this means that the following properties 
hold: 

1. there is a direction in which the flow is strongly contracting ("strongly" means that any 
possible contraction in transverse directions is always strictly weaker); and 

2. transverse to this direction the flow expands areas. 

The robustness of pseudo-hyperbolicity ensures the robustness of the chaotic behavior of orbits 
in the Lorenz attractor. Moreover, pseudo-hyperbolicity is maintained even for (small, indeed) 
time-periodic perturbations so that a periodically forced Lorenz attractor is an example of a gen- 
uinely "strange" attractor since it contains no stable invariant subsets (for more discussion see 
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|SST93j ■ |GUST05] . |TS05j ) . However it is necessary to note that homoclinic tangencies arise in- 
evitably under such perturbations. Consequently systems with periodically forced Lorenz attrac- 
tors fall into Newhouse regions and demonstrate an extremely rich dynamics |GTS93j . |GTS99] . 
Such attractors were called wild-hyperbolic attractors in |ST98j :^ in this paper the theory of 
wild-pseudo-hyperbolic attractor was discussed and an example of a wild-hyperbolic spiral at- 
tractor was constructed. 

Generalizing these ideas, we say that wild-hyperbolic attractors possess the following two 
main properties that distinguish them from other attractors: (1) wild-hyperbolic attractors allow 
homoclinic tangencies (hence, they belong to Newhouse regions); (2) every such attractor and 
all nearby attractors (in the C^-topology with r > 2) have no stable periodic orbits. 

One of the most felicitous properties of wild-hyperbolic attractors is that they can be created 
in local bifurcations of periodic orbits; in particular, when the orbit has three or more Floquet 
multipliers on the unit circle |SST93j . Consequently, the construction of wild-hyperbolic attrac- 
tors in concrete models is possible when the model contains sufficiently many parameters to 
permit such a degeneracy. Indeed, the map is sufficiently complex and has been shown to 
possess wild- hyperbolic attractors |G()STfl5] . Moreover, numerical simulations have shown that 
these attractors persist in quite big parameter domains and can be order-one in size 

In ^we will show that has a small, wild- hyperbolic Lorenz attractor in a (small) param- 
eter domain near the point {Mi, M2, B) = (—1/4,1,1) at which a fixed point with multipliers 
(— 1, —1, -|-1) exists. The general case of fixed points with such a triplet of multipliers was con- 
sidered in |SST93| . where it was shown that two distinct cases occur. One of these cases has a 
local flow normal form that coincides with the Shimitsu-Marioka system (see (|33|) ) . This sys- 
tem has a Lorenz attractor in some parameter domains IShi86j . [Shi93j . As was noted in 
a periodically forced Lorenz system (when the perturbation is rather small) can give rise to 
a wild-hyperbolic attractor. Moreover, this does not destroy the pseudo-hyperbolicity of the 
attractor |ST98j . [GOST05] . Note that the flow normal form only approximates the map and, 
since the residual terms are exponentially small in time, it is reasonable that the map will have 
a wild- hyperbolic attractor (which is small, in fact, because the approximation must be very 
good), see more discussions in {(tOSTOB] . 

The second property of that we mentioned above is that this map can arise as a local 
normal form for multi-dimensional systems with a homoclinic tangency. This shows that such 
maps are not exotic and gives a reason for an extensive study of their properties. 

It is well-known that quadratic maps appear naturally when studying bifurcations of quadratic 
homoclinic tangencies [n??72] . [n?73l| . [T I . Y 8fi| . |( ;TS93j . |( IS' I '03] . They appear as normal forms 
of rescaled first return maps to the neighborhood of a point on the homoclinic orbit. There are 
four standard forms for such one- and two-dimensional "homoclinic maps:" 

1. y = Ml — (the logistic map); 

2. (x, y) = {y, All — Bx — y^) (the standard Henon map); 

3. (x, y) = {y, Ml + M2y - x^) (the Mira map); 

4. [x, y) = (y. Ml — Bx — + £\xy + s^.y'^) (a generalized Henon map). 

^The term "wild" goes back to S. Newhouse |N79| who introduced the term "wild hyperbolic set" for uniformly 
hyperbolic basic sets whose the stable and unstable manifolds (leaves) have a tangency. As Newhouse has proved, 
the latter property is persistent, consequently, systems with wild hyperbolic sets comprise open sets in the space 
of dynamical systems. 
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The third case is a two-dimensional endomorphism and is one of the well-studied polynomial 
maps (see, e.g. |MGBC96] '): we call it the "Mira map" following |GST03j 'l. The fourth case can 
be obtained as a small perturbation (due to e^'^-terms) of the standard Henon map. This map 
was obtained in a study of bifurcations of two-dimensional diffeomorphisms with a homoclinic 
tangency to a saddle fixed point (with unit Jacobian) |GG00j . [GG04j . 

The 3D Henon map is absent from this list, and in ^SJ-^we will fill this gap. 

In ^we study bifurcations of a three-dimensional diffeomorphism with quadratic homoclinic 
tangency to a saddle-focus fixed point. For a three-dimensional map, such points can be of two 
types: (2,1) or (1,2) (sometimes called type A and type B, respectively). Recall that {k,l) 
saddle is a fixed point with dim(T^*) = k, dim(T^") = /, and k + I = n. In our case the two- 
dimensional manifold will be assumed to correspond to a complex conjugate pair of multipliers. 
We will consider the case where the the product the absolute value of the multipliers, Pq is one. 

The bifurcations in the general case, where Pq ^ 1, were studied in jGST96j . |GST03j . We 
will use their "rescaling method" to study our case as well. The essence of this method is that 
a suitable change of coordinates and parameters in a neighborhood of the homoclinic orbit is 
found to bring the map into a standard form. If the homoclinic tangency is assumed to be 
quadratic the resulting map is also typically quadratic. For example, for a type (2, 1) saddle- 
focus with Pq < \, the rescaled first return maps coincide with the standard Henon map up to 
asymptotically small terms (when the time of return tend to infinity). For a saddle- focus of type 
(1,2), the corresponding rescaled first return map is the Mira map. Note the condition Pq < 1 
leads to maps that are effectively two-dimensional; consequently all three-dimensional volumes 
near the fixed point are collapsed. This is not the case when Pq = 1 for the unperturbed map 
when three-dimensional rescaled maps can appear a priori. Indeed, we will show that when 
Pq = 1 the map occurs as the normal form of the rescaled map for type (2, 1) saddle-focus; 
and for the type (1,2) case the corresponding normal form is^ 



x = y, y = z, z = Mi + Bx + M2Z-y^, (2) 

Both and © are quadratic three-dimensional maps with constant Jacobian J = B. It is 
natural to call them the three-dimensional Henon maps of the first and second kind, respectively. 

The quadratic 3D Henon maps are elements of the group of polynomial automorphisms. The 
study of the dynamics of polynomial mappings has a long history applied dynamics; for example, 
they often used in the study of particle accelerators |DA96j . The study of Cremona maps, i.e.. 



polynomial maps with constant Jacobian, also has intrinsic mathematical interest Eng58 . An 
interesting mathematical problem is to obtain a normal form for an arbitrary degree-d Cremona 
map. This problem is unsolved; the major obstruction is that it is not known in general that such 
a map has a polynomial inverse. This is the content of the much studied, "Jacobian conjecture": 

Conjecture (O.T. Keller (1939)). Let / : ^ M" be a Cremona map. Then f is bijective 
and has a polynomial inverse. 

This conjecture is still open. 

The set of polynomial maps with polynomial inverses is called the affine Cremona group. 
For the case of the plane, the structure of this group was obtained by Jung |Jun42j . Jung's 
theorem states that planar polynomial diffeomorphisms are "tame;" that is, they can be written 

^Note that ij^Jl is also the rescaled first return map for a four-dimensional diffeomorphism that has a quadratic 
homoclinic tangency to a saddle-focus fixed point with two pairs of complex conjugate leading multipliers 
l(jT«93l . inHT96l . kiHT03l . 
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as a finite composition of "elementary" (or triangular) and affine maps. More recently Friedland 
and Milnor |FM89j showed that any map in this group is either conjugate to a composition of 
generalized Henon maps, 

X = y , y = -Bx+p{y) , 

where p{y) is a polynomial and B ^ 0, or else the map is dynamically trivial. 

Not much is known about the structure of the affine Cremona group for higher dimensions, 
though it was recently proved that the three-dimensional Cremona map 

X = X — 2y{xz + y^) — z{xz + y^)^ , 
y = y + z(xz + y2) , 
z = z , 

is not tame |SU04j . Thus the normal form for the 3D affine Cremona group must consist of more 
than tame maps. However, it is still possible that the nontame maps are dynamically trivial, 
like the above example. 

Another interesting problem concerns the approximation of smooth diffeomorphisms by poly- 
nomial diffeomorphisms. One result along this direction is for the class of C°° symplectic diffeo- 
morphisms, which, as Turaev has shown |Tur03j . can be approximated by compositions of the 
symplectic Cremona maps 

X = Xo + y , y = -x + Vp{y) , 

where x, y, xq £ M" and p is a polynomial. Related approximation methods are used in the study 
of particle accelerators |DA96j . and results have also been obtained for analytic symplectic maps 

Some work has been done on finding normal forms for quadratic symplectic and quadratic 
volume-preserving maps. Moser showed that every quadratic symplectic map in M^" can be 
written as the composition of an affine symplectic map and a symplectic shear |Mos94j . For 
the four-dimensional case, he obtained a normal form. For the case of volume-preserving maps, 
it was shown in |LM98j . |LLM99j that any quadratic three-dimensional, volume-preserving map 
whose inverse is also quadratic and whose dynamics is non-trivial can be transformed (by affine 
transformations) to the form 

X = y, y = z, 

(3) 

z = Mi + Bx + Mly + M|z + ay^ + hyz + cz^ , 

where B = 1, |a + 6 + c| = 1 and one of the parameters or Ml can be eliminated (by a 
coordinate shift). By analogy with the planar results, there are two additional cases but they 
have trivial dynamics (see We also note that it is not hard to construct three-dimensional 
quadratic diffeomorphisms that do not fall into this classification because they have an inverse 
of higher degree. For example, the diffeomorphism 

x = x, y = y + x^, z = z + y'^ , 

has a quartic inverse. 

Note that both of the maps ((TJ) and © fall into the form ©. Thus the results of |LM98j 
show that when B = 1 both of the 3D Henon maps occur as normal forms for quadratic elements 
of the affine Cremona group. Moreover, even when B ^ 1 the map Q is a normal form for a 
certain class of quadratic maps [TatOlj . Tatjer showed that every quadratic map, /, for which 
/ o / is also quadratic, but has no invariant linear foliation is linearly conjugate to 0. In ^ 
we will generalize the results of |LM98j to show that Q is the the normal form for quadratic 
diffeomorphisms with quadratic inverses when B ^ 1 as well. 
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2 Statement of Results 

Consider a three-dimensional diffeomorphism go £ C(M^) for some r > 5 that satisfies the 
fohowing conditions (see Fig. : 

(A) gQ has a type (2, 1) saddle-focus fixed point O with multipliers (Aoe*"^", Aoe~*'''°,7o) where 
< Ao < 1 < ItoI (w.l.o.g. < < tt),; 



(C) The stable, W^, and unstable, W^, invariant manifolds of the fixed point O have a 
quadratic homoclinic tangency at the points of a homoclinic orbit Tq. 



Figure 1: Geometry of the homoclinic orbit to a type (2, 1) saddle-focus for the map go under assumputions 



Diffeomorphisms that are close to go and satisfy conditions (B) and (C) comprise, in the 
space of diffeomorphisms, a codimension-two, locally-connected surface Tl — the bifurcation 
surface. 

Our main problem is to study bifurcations of periodic orbits in families of diffeomorphisms 
that transversely cross TC. We denote such a family by ge so that go satisfies (A)-(C). The first 
question that arises is the optimal choice of the governing parameters e. It is evident that the set 
of such parameters must include both a parameter fi = £i of the splitting of the initial quadratic 
homoclinic tangency and a parameter £2, similar to £2 = Pg — 1, that varies the Jacobian of the 
saddle-focus fixed point. However, it was shown in [(TonODj that diffeomorphisms on Tl possess 
continuous (topological) conjugacy invariants on the set of non-wandering orbits, so-called 
moduli. The most important modulus is the angular argument ip of the complex multipliers of 
the saddle-focus (see By definition, O-moduli are natural governing parameters. Therefore, 
we include the S7-modulus if into the set of governing parameters. Thus, we take as the third 
parameter £3 = (p — ipo- 

More generally, we must consider a three-parameter family g^, e = (£1,62, £3), that is 
transversal to at £ = and such that 



Consider a sufficiently small and fixed neighbourhood U of the orbit Vq. It can be represented 
as a union of a small neighbourhood Uo of O with a number of small neighbourhoods of those 



(B) Po 



Jol = IA270I = 1. 




W"(0) 



A 



(A)-(C). 
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points of Fq which do not belong to Uq. In the family we study bifurcations of single-round 
periodic orbits from U — that is orbits periodic that enter Uq exactly once. Every point of such 
an orbit of period k can be considered as a fixed point of the corresponding first return map 
g^ along orbits lying in U. In ^ we will prove the following result about the form of this first 
return map in some domains of the parameters. 

Theorem 1. Let g^ be the three parameter family defined above. Then, in any neighborhood of 
the origin, e = 0, of the space of parameters, there exist infinitely many open domains, A^, such 
that the first return map g^ takes the form 

X = y + o(i)fe^+oo , 

y = z + o(i)fe^+oo , (4) 

Z = Mi + BX + M2Y - Z2 + o(l)fe^+oo , 

where {X, Y, Z) are rescaled coordinates and the domains of definition of the new coordinates 
and the parameters {Mi, M2, B) are arbitrarily large and cover, as k ^ 00, all finite values. 
Furthermore, the Jacobian, B, is positive or negative depending on the orientability of the map 
g^. The terms denoted o(l) are functions (of the coordinates and parameters) that tend to zero 
as k ^ +CXD along with all derivatives up to order r — 2. 

Theorem^shows that the 3D Henon map (pQ) is a normal form for this homoclinic bifurcation, 
since it appears when we omit the o(l)-terms in Q). We will see that the new parameters Mi, M2 
and B can be given as functions of e = {fi,ip — fo, Ps — 1) by 

Ml ~ 73'=(/x + 0(A^)) , M2 ~ (At)'^ cos(A:^ + /?) , ^ ~ (A't)'' 

where P is some coefficient (an invariant of go) and (A(0),7(0)) = (Ao,7o)- The exact relations 
will be given in formulas pflj) . (|25|) and (|27)) below. Thus, we can study the bifurcations of single 
round periodic orbits in the family g^, when e E A^,, by studying the bifurcations of fixed points 
of the 3D Henon map (^. 

Certainly, the problem of studying bifurcations in the 3D Henon map is much simpler than 
the corresponding problem for the original first return maps. Moreover, the 3D Henon map 
provides a simple, standard model for which various well-tried methods (including numerical 
ones) can be applied. However, the bifurcation problem as whole for the 3D Henon map is very 
complicated and here we focus only on two aspects. In ^we first give the simple construction 
of bifurcation surfaces, lines and points corresponding, respectively, to codimension one, two 
and three bifurcations of the fixed points. Then we consider the problem of the existence of 
wild-hyperbolic strange attractors. We will obtain the following result. 

Theorem 2. Let g^ be the three parameter general family defined as before. Then, in any 
neighborhood of the origin e = of the space of parameters, there exist infinitely many open 
domains 6^ C A^. such that the diffeomorphism g^ of the family at e £ 6k has a wild-hyperbolic 
strange attractor of Lorenz-type. 

The proof of this theorem is given in gni following the analysis of |(;()STn5j . The main idea 
of the proof is to apply the results of paper |SST93j where it was shown that maps having 
a fixed point with the triplet (—1,-1,-1-1) of multipliers can, in some cases, be reduced to a 
flow normal form that coincides with the Shimitsu-Marioka system (see (|3f)|) ) which possesses a 
Lorenz attractor for some parameter domains |Shi86j . [Shi93j . Our verification of these facts for 
the corresponding fixed point of the 3D Henon map allows us to obtain Th. [3 A justification 



7 



of the fact that the attractor in the 3D Henon map is actually the wild-hyperbolic Lorenz-type 
attractor is given in |GUST05j . Consequently, we will omit these details. 

Solving the bifurcation problem for the family allows us to solve almost automatically 
an analogous problem for a diffeomorphism G C (r > 5) with a type (1,2) saddle focus. 
Specifically, assume that qq satisfies (see Fig. 

(A') go has a type (1,2) fixed point O with multipliers {Xo,joe^^°,'joe~^^°) where < |Ao| < 
1 < 70 (w.l.o.g. < -00 < vr); 

(B') Po = |Ao7ol = 1; and 

(C) the stable and unstable invariant manifolds of O have a quadratic homoclinic tangency at 
the points of some homoclinic orbit Tq. 




Figure 2: Geometry of the homoclinic orbit to a type (1,2) saddle- focus for the map go under assumptions 
(A')-(C'). 

As in the case of the type (2, 1) saddle-focus, we consider a three-parameter family g^ where 
e = (/i, — 1, ■0 — ipo)- We suppose here that /i is a splitting parameter. Pi = |A(e)7(e)^| and 
A(0) = Ao,7(0) =70- 

Theorem 3. Let g^ he the three- parameter general family defined above. Then, in any neigh- 
borhood of e = there exist infinitely many open domains, A^, such that the first return map 
g^ takes the form 

X = Y + o{l)k^+^, 

y = Z + o(l)fc^+oo, (5) 
Z = M, + BX + M2Z - y2 + o(l)fe^+oo, 

where (X, Y, Z) are rescaled coordinates and the domains of definition of the new coordinates 
and the parameters {Mi, Ad2, B) are arbitrarily large and cover, as k ^ 00, all finite values. 
Furthermore, the Jacohian, B, is positive or negative depending on the orientability of the map 
g^.The terms denoted o(l) are functions (of the coordinates and parameters) that tend to zero 
as k ^ +00 along with all derivatives up to order r — 2. 

Proof. The proof immediate from Th.|2 Indeed, the inverse gQ^ is, evidently, a diffeomorphism 
satisfying the original conditions (A-(C). Thus, the map ^ is, in this case, the rescaled map 
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for the inverse first return map g- . The last can be be easily written in form @, as we show 
explicitly below. □ 

The inverse of is of interest in any case. Since has the constant Jacobian B it is 
invertible whenever B ^ 0. This inverse is 

y = x, z = y , 

x = ^{-Mi + z-M2y + y^) . 

This map can be rewritten in the "standard" form of |LM98j . if we define new variables 

^^new ^yuew ^ ^new^ _ B~^{z,y,x) and new parameters 

then the resulting map has precisely the form @. As we stated above map is well-known in 
homoclinic dynamics: it was found to be the rescaled first return map for four-dimensional 
diffeomorphisms that have a quadratic homoclinic tangency to a saddle-focus fixed point with 
two pairs of complex conjugate leading multipliers |GTS93j . |GST96j . |GST03j . 

The dynamics of each of the maps and is extremely rich. However, for the volume- 
contracting case, we note that these maps demonstrate distinct types of stable dynamics. For 
example, gives rise to attractors similar to those of the Lorenz system, as we will see in 
^ However, need not have attractors of this "wild" type, although certainly, a "Lorenz 
repeller" will be present. 

Nevertheless, we can expect existence of "extremely small" Lorenz-type attractors in Q by 
the following argument. As we will see in ^ the map ((2j can have a fixed point with multipliers 
(-|-1, e*"^, e~*'^) at certain values of its parameters. The fiow normal form near this point can 
have a homoclinic loop to a saddle-focus and, thus, ((2)) can have a fixed point O with multipliers 
(7, Xe'^, Xe-'^), where < A < 1< 7 and = 1, such that the manifolds W'{0) and W^iO) 
will have a (quadratic) homoclinic tangency. Consequently, we can apply Th.l^Jto this situation 
giving immediately the promised "extremely small" attractors. 

Nevertheless, as it is shown in |GUST05] the Lorenz-type wild-hyperbolic attractor for 
can be quite big and were observed to have a large basin. The proof of existence of big, globally 
attracting wild-hyperbolic attractors in map Q is still open problem. 

Finally, we consider the more general structure of polynomial diffeomorphisms of M^. As we 
remarked in the introduction, the 3D Henon maps and © have the form of Q when B = \. 
As we will shown in ^ this result has the following generalization; 

Theorem 4. Let / : M'^ — > 6e a quadratic dijjeomorphism with a quadratic inverse and 
constant Jacobian det Df{x,y, z) = B. Then, after an affine coordinate transformation f can 
be written in the form ^ or else the dynamics of f is trivial. 

Note that both of the maps Q and ((2)) fall into this classification. The classification of the 
dynamics of the trivial cases can be done from the normal forms H43() that will be obtained in 

m 

3 Local and global maps 

We begin the proof of Th.^by constructing a local map Tq, valid in a neighborhood Uq of the 
fixed point O, that simplifies the description of the local unstable VF;"^(0) = W^{0) n Uq and 
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local stable Wi^^^{0) = W'^{0) D U manifolds of O. Then we obtain a global map Ti that takes 
a neighborhood of a point on Wl^^{0) to a neighborhood of a point of Wf^^{0). 

By |GS93j . [5STC98j . one can introduce coordinates {x,y) = {xi,X2,y) in a neighbor- 

hood U of the fixed point O such that the local map Tq^e) = g^^Ug can be written, for small 
enough e, in the form^ 

X = \R^x + i{x,y) 
^0 : _ . , ^ (6) 



where 



y = 7y + r]{x,y) 



\ sin If cos (/? / 



is the rotation by angle f, and 

\m=0{\\xfy), andr] = 0{\\x\\y^) . (8) 

Here the coefficients, A, 93, and 7 as well as ^ and 77 depend (smoothly) on the parameters e. 
The form Q will be called the main normal form; it is very convenient for calculations. In 
particular, the local stable and unstable manifolds of O are Wf^^ = {{x,0)} and W^^^ = {(0, y)} 
in the neighborhood U. 

When e = 0, choose a homoclinic point = (0,0, y^) G W^^^ on the homoclinic orbit Tq 
and a corresponding point = {xf,X2,0) G ^foc that is an image of M~: = gQ°{M~) 
for some uq > 0, recall Fig. ^ Let 11^ and 11" be some small neighborhoods of the points 
and M~ . The global map Ti{e) = Qs'^^yi- '■ ~^ defined for small enough e, and 

can be written in the coordinates of © as 

^^ x-x+ = Ax + b{y-y-) + 0{\\xf + \y-y-\'^), 

y = y+(^s) + c^x + d{y-y~)'^ + 0{\\x\\^ + \\x\\\y-y-\ + \y-y~\^), 

where the 2x2 matrix A, the vectors x"*", b, and c, the coefficients y"*" and d, and the higher order 
terms all depend smoothly on e. We noted especially the dependence of the coefficient y'^is) 
that will become the splitting parameter for the manifolds and near the homoclinic 
point M+. As discussed in the previous section, we therefore denote it by 

^ = y+(e) (10) 

and consider as the first of the governing parameters. 

The assumption (C) that the tangency between and is quadratic implies that 

dy^O. (11) 

In addition, since go is assumed to be a diffeomorphism, and since the Jacobian 




DT,{M-) ={ ^ J , (12) 

then we must have ||5|| 7^ 0. 



^It follows from pS93 , SSTC98J that © holds for general case, including Pq = 1. For the volume-preserving 
case the analogous form was established in IBaz93| 
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Indeed if 62 (^) is not identically zero, then one can use the linear coordinate change, 

where ii^, 0, is a rotation by an angle to transform the map Q to a new map of the 

same form such that the vector b becomes 

6^ = (6i,0). (14) 

Here the new matrix point x"*", and vector c are modified from the corresponding coefficients 
in the original ©• 

Now the nonvanishing of the determinant of the Jacobian (|12)) implies that 



Ji = detL'ri(M-) = 6i(a2iC2 - 02201) / . (15) 

One of merits of the main normal form ^ is that the iterations Tq : [/q — > Uq for any k can 
be calculated quite easily, namely, in a form close to the form when Tq is linear. In particular 
for the linear case we have 

Xk = X^Rk^XQ , yo = I'^Vk ■ 

This is the so-called cross-form of map T^. We can write Tq in an analogous cross-form when 
Tq is in the form (jB)) 

Lemma 5. For any positive integer k and for any sufficiently small e the map Tq (e) : (xq, yo) — > 
ixk,yk) can be written in the cross-form 

Xk = X'^Rk^xo + X^''Ck{xo,yk, e) , 

_. , _, (16) 

2/0 = 7 Vk + Xj r]k{xo,yk,e) , 

where the vector S^k '^^'d the function rjk are uniformly bounded along with all derivatives up to 
order (r — 2). In addition, all of the derivatives of order (r — 1) for the functions Xk and yo tend 
to as k ^ 00 . 

Proof. Let {xi,yi) £ U , i = 0,...,k, denote the orbit: {xi,yi) = To(xj_i, yj_i). Then it is easy 
to see that the j^^ iterate of (jH)) is 

Xj = X^Rj^xo + ^ A-'"''"^i?(j_s-i)^?(2;s,ys,e) , 



yj = 7^ ''yk-J2^^ ' ^r]{xs,ys,e) 



S=] 

The solution of ((T7j) for any given sufficiently small xq, yk and large k (the so-called boundary 
value problem) can be found by successive approximations |AS73j if the zeroth approximation 
is taken as the corresponding solution of the linear problem. Using the fact that IA7I > 1, we 
obtain the following estimates 

\\xk - X^Rk^xoW < L e'a'^-^-U^^ItI^-^ < TA2^ 

s=0 

iyo-7'Vi <i i;'i7r'"'Ai7i''-'' <^A'=i7r', 

s=0 
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where L and L are some positive constants estimating the norms for derivatives of the functions 
^ and ry in the domain {/, from (jH)). These inequahties evidently imply the desired result for 
the map. Estimates for the derivatives can be found analogously. General results of this form 
for near saddle fixed point were given in [GS93j , ^STC9&| . Here we give sharper estimates that 
apply in this particular case. □ 



4 Rescaled first return map: Proof of Th. [T] 

The main task of this section is to obtain an analytical form for the first return map for diffeo- 
morphisms of the family to a neighborhood of the homoclinic point M+. To construct this 
map we first find an image of the local map Tq for which Tq (11^) nil" 7^ 0. This will be true for 
each k >k where k is some sufficiently large integer. For each such k there are strips cr^ G 11+ 
and (7^ C n~ such that ro((T^) = <7^. A first return map is obtained by composing Tq with 
the global map Ti : 11+, to obtain a map on a neighborhood of M"*", recall Fig. Q Thus 

the first return map will be of the form 

T' rp rph 

k = ^l-tQ ' 

Under the global map Ti the strip is transformed into three-dimensional horse-shaped region 
Ti(cr^) that has non-empty intersection with when e is small enough. Thus, the map is 
a generalized, three-dimensional, horseshoe map. We will show that these horseshoe maps can 
possess rich dynamics; for example, they can Lorenz-like, wild-hyperbolic strange attractors. 
The goal is the prove the following lemma: 

Lemma 6 (Rescaling Lemma.). For any sufficiently large k such that 

A21 (kip) = 021 cos kip + 022 sin kip ^ (18) 

the first return map T}^ can he brought, by means of an affine transformation of coordinates, to 
the following form 

Y = Z + 0{\^) , 

X = Y + 0{\^ + \H'^), 

Z = Ml- Z'^ + hi\-i^{ci COS kip + C2s\nkip + vl)Y 
+61^21 A2^7'=(c2 cos kip - ci sin kip + vl)X + 0{\^) , 

where v^'^ = 0(A'^) are some coefficients, 

Ml = -d72^[/i - 'j-'^y- + X'' {coskip{cixfc2X+) + smkip{c2X+ - C2X+)) + 0{X^'')]. (20) 

Also, in my\) the O-terms have the indicated asymptotic behavior for all derivatives up to order 
(r-2). 

Proof. Using formulas Q, (|16|) and conditions (|1()|) and (|14j) we can write the first return map 
Tfc = TiTq in the form 

Xi 
X2 

-f-^y + \''-i-^r}k{x,y,e) 



= xj + \'^{eiiXi+ei2X2) + hi{y-y ) + ai, 

= X^ + \''{e2lXi+ e22X2) + OL2, ^^^^ 

= ^ + \^{ci cos kip + C2sm.kip)xi + 

X^(c2 cos kip — ci sin kip)x2 + diy — y~)^ + 03, 
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where the 2x2 matrix E = ARi^^p and 

ai,2 = C'(A2'=||x|| + A'^llxllly - y-\ + {y - y'f), 
as = ©(A^^-'llxll + A'^llxllly — y^\ + \y — 

Introduce new coordinates 

where Pk = C>{X^), in order to ehminate the affine terms in the first two equations and the term 
hnear in y in the third equation of (|21|) . Then, the map H21() is rewritten as 

xi = eiiX'^xi + ei2X''x2 + by + O {X^'' \\x\\ + A^'llxllllyH + ||y|p), 
X2 = €21X^x1 + €22X^x2 + ©(A^^'llxll + A''||x||||y|| + ||y|p), 

y + A'=0(||x|| + ||y||) = M + d7V+ (22) 
+A^7^(ci cos kip + C2 sinkip + 0{X^))xi + X^^^{c2 cos kip — ci svukip + 0{X^))x2+ 
+-f^O{X^^\\xf + X^\\x\\\\y\\ + ||yf ), 

where dnew = d + 0{X^) and 

M = -^''[^ - j~''y~ + X'^ [coskip{cixfc2X2) + smkp{c2xf - C2X^)) + 0{X'^'^)]. (23) 

Next, we rescale the map introducing new coordinates (X,Y,Z) defined through 

-1 = -^il-'Y, X2 = -^AS-^'X, y = -W^Z, (24) 
d d d 

where we suppose that values of 99 are such that ()18|) holds. Then, map ()22() is rewritten in the 
new coordinates in form (|19|) . □ 

Now theorem ^ easily follows from this lemma. Indeed, we have only to indicate those 
domains of the initial parameters where the coefficients in the third equation of (|19j) are 
finite. Denote the coefficient of the third equation of ()19() in front of Y as M2, i.e. 

M2 = hi X^^^ (ci cos kip + C2sm.kLp + v\). (25) 

Since I7I > 1 and IA7I > 1, the coefficients Mi and M2 can take arbitrary finite values when 
varying parameters // (for Mi) and ip (for M2), see formulas (|20|) and H25() . Moreover, the 
finiteness of M2 means that the corresponding values of the trigonometrical term ci cos kp + 
C2sin/c(/9 are asymptotically small and, in any case, do not exceed |A7|~^ in order. For such 
values of we obtain 



C2 cos kp — ci sinkp = iby^cf +03(1 + ...), 

A - 7 , -I, ,a2lC2-a22Cl (26) 

^21 021 COS A:(/3 + 022 sm fcyj = zb , (1 + ...J ^ ^ 

Thus, by (Cni), A21 7^ in this case. Denote as B the coefficient in front of X2 from the third 
equation of ()19() . By virtue of H26|) and (jlSf) we have 

5 = 6i(a2iC2 - a22Ci)A2V(l + •••) = JiA2V(1 + •••)• (27) 

Since |A^7| = 1 at e = 0, it follows from (|27|) that -B can take (when A; is large) any arbitrary 
finite value (positive or negative depending on the sign of Ji}?^^^) as the value of the Jacobian 
A^7 of the local map Tq varies. This completes the proof of theorem ^ 
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Remarks 



1. We noted that at least three parameters must be used to study bifurcations of a homochnic 
orbit for a three-dimensional saddle- focus. The governing parameters that we use are 
ei = n corresponding to the splitting distance between the stable and unstable manifolds, 
the phase £2 = ~ V'o of the focus multiplier, and the Jacobian £3 = |A^7| — 1 of the 
saddle- focus. 

2. If the map go is exactly volume preserving, i.e., det Dqq = 1, then the resulting map ((T^ is 
as well, i.e. B = 1. Such maps were studied in |LM98j . Moreover, Th.|21implies that there 
are nearby maps (by unfolding the parameter B to make the map volume contracting) 
that have strange attractors. 

5 Local Bifurcations of 3D Henon maps 

In this section we will study the dynamics of the 3D Henon map (pQ) and prove Th. |2 The 
dynamics can be complex, especially when |-B| is close to one, essentially because when |S| = 1 
the map has, for certain values of the parameters, fixed points with three multipliers are on the 
unit circle. It is well known that bifurcations of such points can lead to appearance of homoclinic 
tangencies, invariant circles, and even strange attractors (as \B\ decreases below one) or repellers 
(as \B\ increases above one) |Shi81j . |?^ST93j . 

We are interested primarily in the stable dynamics of the map (^, i.e., in characterizing its 
attractors. One way to find attractors with complex dynamics is to study the local normal forms 
that arise near fixed points that have three multipliers that are ±1. We follow here the analysis 
of |SST93j where local normal forms of the equilibrium of a flow with three zero eigenvalues and 
with additional symmetries were analyzed and conditions for the existence of strange attractors 
were given. ^ One of the most interesting of their results is that Lorenz attractors can appear 
by local bifurcations of an equilibrium with three zero eigenvalues when there is a reflection 
symmetry through the origin |Shi81j . [5ST93j . [PST98| . These results can be applied to the local 
bifurcation analysis of when the fixed points has multipliers (+1, —1, — !)• 

As is true for the general volume-preserving quadratic diffeomorphism of Lomeli and Meiss 
|LM98j . the 3D Henon map has at most two fixed points on the diagonal: 

p± = (4.4.4) . 

1 / ^\ (28) 

z± = -[-l + B + M2±VAj , 

providing that the discriminant, A, is nonnegative, 

A = (M2 + S - 1)2 + 4Mi > . (29) 

The primary bifurcations of these fixed points occur when they have multipliers on the 
unit circle. These bifurcations are easily determined using the characteristic polynomial of the 
Jacobian of 

p{X) = + 2zX^ - M2A - B . 

^These results can be applied to the analysis of bifurcations of codimension-three fixed points of three- 
dimensional maps. In this case we can often construct some approximated fiow normal form whose time-one shift 
map coincides (up to terms of a certain order) with some power of the original map restricted on a neighborhood 
of the fixed point. Usually, the presence of multipliers —1 or e^"^ implies that the corresponding approximated 
flow possesses some symmetry, such as a reflection x — * ~x, or a rotation, 02- 



14 



The first codimension-one case corresponds to a single multiplier one, where p{l) = 2z* + 1 — 
B — M2 = 0, or equivalently to the surface where A vanishes: 

L"" : Ml = -^{M2 + B -if . (30) 



This is simply a saddle-node bifurcation that creates the pair of fixed points (|28|) in the domain 
A > 0. Similarly, a period-doubling bifurcation occurs when there is a multiplier —1, i.e., where 
p{—l) = 2z* — 1 — B + M2 = 0, or equivalently for 1 — M2 = ±|\/A, which gives 

LP'^ : Mi = ^{M2-l- B){3M2-3 + B) . (31) 

Note that this bifurcation occurs for the point P+ when M2 < 1, and for the point P_ when 
M2 > 1. 

The final codimension-one bifurcation corresponds to the Hopf bifurcation when pair of 
multipliers is on the unit circle, A = e^**^, so that 

p{X) = - (5 + 2 cos ip)X^ + {l + 2B cos ip)\-B , 

or equivalently on the surface 

. i Mi = {cos^ + lB){{l-2B)cos^ + lB-2), 
' \ M2 = -l-2Bcos(p . 

Here we have written the equation for the surface L"^ in parametric form using the phase < 
ip < IT. When B < 1, this bifurcation only occurs for the fixed point P+. 

There are three codimension-two curves in the {Mi, M2, B) parameter space. The Hopf 
surface, terminates at L^" when the multipliers are (-6,1,1) along the curve (Mi,M2) = 
{j{B + 2)^, —2B — 1). It also terminates on L^"^ where the multipliers are {B, —1, —1) along the 
curve {Mi,M2) = {\{B — 2){7B — 6),2B — 1). Finally, the period-doubling and saddle-node 
curves are tangent when the multipliers are {—B, —1, 1) on the curve (Mi, M2) = {—^B^, 1). 

The bifurcation surfaces are most easily viewed in two-dimensional slices through constant 
B planes. An example of the resulting bifurcation curves are shown when B = 0.6 in Fig. |31and 
Fig. 01 Note that P+ is stable inside the curvilinear triangle bounded by L'^", LP'^ and L'^, when 
B < I. The second fixed point, P_, is never stable when B < 1. 

A similar analysis can be used to derive the bifurcation surfaces for the map (jSJ. Indeed, 
the formula for the fixed points, and consequently the saddle-node surface is the same as for 
The period-doubling and Hopf surfaces are slightly modified. An example of a slice through 
these surfaces for B = 0.6 is shown in Fig. Eland Fig. (HI 



6 Wild-hyperbolic Lorenz attractors: Proof of Th. [2] 

We will now focus on the point P+ and expand around the codimension-three bifurcation point 
where the multipliers are (—1,-1,1), that is about the point {Mi, M2, B) = {—^,1,1) where 
z* = 1/2. Note that on the parameter slice B = 1, the Hopf curve H32|) coincides with the 
saddle-node curve H30() so that there is no Hopf bifurcation. Moreover the codimension-two 
points in Fig. |31 corresponding to the multipliers (— 1,— 1,-B) (square) and (— l,l,i?) (triangle) 
collide at the codimension-three point. Otherwise the bifurcation diagram for P+ looks similar 
the slice in the figure. 

Our goal is to prove Th. (21 This can be done with the aid of the following lemma. 
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Figure 3: Bifurcation curves of the map Q for the fixed point P+ when B — 0.6. Fixed points do not exist in the 
shaded region. The solid curve is the saddle-node L"", the dotted curve is the period doubling, L*"*, and dot-dash 
curve is the Hopf bifurcation, L'^ . The fourth curve (dashed) corresponds to the existence of a double eigenvalue. 
There are five corresponding regions, and a representative configuration of the multipliers in the complex- A-plane 
is shown for each. The codimension-two bifurcations are shown as the "square" (period-doubling- Hopf ) , "triangle" 
(saddle-node-period doubling) and "'circle" (saddle-node-Hopf) points. 



Lemma 7. In a region of the parameters (Mi, M2, B) near the point (—1/4,1,1) the second 
iterate of in some neighborhood of -Pf is approximately the time-one map of the the three- 
dimensional flow 

x = y , 

y = x{l - z) - Xy , (33) 
z = —az -\- , 

up to arbitrary small correction terms. Moreover the coordinates (x, y, z) and parameters A and 
a can take arbitrary finite values. 

Proof. Shifting coordinates x ^ x — Zj^,y ^ y — Zj^^ z ^ z — Zj^ brings map to the form 

X = y, y = z, z = Bx -\- M2y — Izj^z — z^ (34) 
Introduce new (small) parameters ei,e2)e3 by the following formulas 

ei = l-B, £2 = 1- M2, £3 = 22+ - 1. 

Note that 

ei > 0, ei+e2 + e^ = VA > 0. 
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-3.0 0.0 3.0 

M, 

Figure 4: Bifurcation curves of the map Q in for the fixed point P_ when B = 0.6. The curves have the same 
meaning as in Fig. |21 For this fixed point there is no Hopf bifurcation. 

The first inequality means that B < 1, and the second is fulfilled automatically when the fixed 
points exists, e.g. to the right of the L^" curve in Fig. |21 We have also the following relation 

l+4Mi =ei + 2(l + e3)(ei+e2) 

between the new parameters and the parameter Mi. Next, we bring map (|34j) to the following 
form 

Ul = —Ul + U2, 

3 -i 

U2 = aiui + (-1 + a2)u2 + U1U3 - -U2U3 + 0(11^11 ), (35) 
11 

% = (1 + a3)u3 - -ul - -ul + Oi\\uf), 

where 

2qi = £1- 82 + £3, 4a2 = ei + £2- Sea, 4a3 = -{ei + £2 + £3)- 

Note that the Jacobian of the map H35|) is close to the Jordan normal form for the multipliers 
(— 1, —1, 1) when Qi, 03 are small. The first equation of 1)35(1 does not contain any high order 
terms, and in the second and third components, we take into account only quadratic terms. 
Moreover, these quadratic terms have been further simplified by standard normal form analysis, 
so that only the resonant ones remain. 

The second power of the map ()35|) can be embedded into a flow. Namely, using the Picard 
iteration procedure (of the second order), we can see that the map (|!^ coincides, up to terms 
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Figure 5: Bifurcation curves of the map l|5J for the fixed point P+ when B = 0.6. The curves have the same 
meaning as in Fig. |21 



of order Odltip), with the time-one map of the following flow 

5 

ill = -2^2 - UiU^i + 





U2 = PlUi - /32U2 - 2uiU3 + 2U2U3 (36) 

1 2 1 2 1 2 1 

U3 = P3U3 - -Ui - —U2 - -U3 - -UlU2. 



After a smooth change of coordinates of the form 

vi = ui, V2 = iti, f 3 = ns + Odlnll^) 
equations are written as follows 



where 



Vi = V2 

V2 = P-vi - Av2 - aviV3 - aiV2V3 (37) 
V3 = -a + vl + b{vf + vj), 



3 „ , 1 _ „ ;r ei + 3e2 - 3^3 _ (£1+^2 + £3)^ 
ai = 2, 0= -, fi = -El - 362 - £3, A = , a = — 
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Figure 6: Bifurcation curves of the map in for the fixed point P_ when B = 0.6. The curves have the same 
meaning as in Fig. |21 For this fixed point there is no Hopf bifurcation. 



Flow (|57j) is the same that was considered |SST93j and, in our case we have also that ab > 0. 
Then, we can introduce new time r and new rescaled coordinates {x,y,z), by formulas 



r = st, where = /i + aVoi = ei — £2 + 



£3, 



(38) 



Vl 



X\l — , V2 

ab 



V3 



-\fa H 

o 



such that system (|57j) is reduced to the following one 



(39) 



X = 7/, 

y = x(l-z) -Ay + 0(s), 
i = -az + + 0(5), 



where 



a 



£1 + £2 + £3 



£1 + 3£2 - 3£3 



(40) 



(41) 



r 2^£i - £2 + £3 ' T 2^£i - £2 + £3 ' 

Since values of s can be arbitrarily small, it implies that the rescaled coordinates (x, y, 2) and 
new parameters A and a can take arbitrarily finite values. □ 

Note that if we omit the ©(s) terms in the system ()4U() . we obtain system ()33|) . which is 
the well-known, Shimizu-Morioka model. This system was intensively studied in jShi86j . [Shi93j 
where, in particular, the existence of the Lorenz attractor was proved for some domain of the 
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parameters. For example, the domain of the existence of the Lorenz attractor in Shimizu- 
Morioka model is known to be contained in a rectangle 0.7 < A < 1.5, 0.1 < a < 0.6 f |SST93j . 
fig. 13). Using this fact we can pointed out the following boundaries of values of the parameters 
(ei, £2, £3): 

£1 > , 

£1 + £2 + £3 > , 
£i-£2 + £3>0, 

£l+3£2-3£3 ^ ^ . ^42) 

0.7 < — I < 1.5 , 

2V£i -£2 + £3 

0.1 < ^ ;+"^ + ^^ < 0.6 , 
2V£i - £2 +£3 

(and, respectively, of the parameters (Afi, M2, .B)) such that in the corresponding domain of the 
initial parameters (Mi, M2, -B) our 3D Henon map has a wild-hyperbolic Lorenz-type attractor. 
This completes the proof of the main theorem. 



7 Quadratic Volume-Preserving Mappings: Proof of Th. HI 

Suppose that / : ^ is a quadratic Cremona map, i.e., is a diffeomorphism with constant 
Jacobian det(D/) = B. Then it is easy to see that without loss of generality we can write it in 
the form 

/ = 5 o L o (id + Q) . 

where S{x^y^z) = (B^^^x, B^^^y, B^^^z) is the diagonal scaling map with det(S') = B, L is an 
affine volume-preserving map, Q is a vector of quadratic forms so that the map x + Q{x) is a 
volume-preserving quadratic shear |LM98j . 

As in Th.|l]we now assume that the map / has a polynomial inverse (i.e., it is an element 
of the affine Cremona group) and that the inverse is also a quadratic map. The same property 
then will also hold for the volume-preserving map g = o f. It was shown in jLM98j that 
maps of the form g can be transformed, by an affine change of variables, to one of the three 
following forms: 

gi{x,y,z) = {y,z,Mi + M2Z + M^y + x + q(y, z)) , 
g2{x,y,z) = {xo + Jf^^^Vo - Miz,zo + M2Z + y + q{x, z)) , 

93{x,y,z) = {xo + y + Mix,yo- -^x, zo + M2Z + q{x,y)) , 

IVI2 

where q{x, y) = ax^ + hxy + cy"^ is a scalar quadratic form. 

Consequently g = A o g-i o for some affine transformation A, and i = 1, 2, or 3. Using 
the same transformation on / gives 

f = A-^ofoA = {A-^ oSoA)ogi. 

But the scaling map S commutes with A, so 

f = S ogi . 
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Finally we can apply an additional scaling transformation to the variables and parameters to 
transform / to one of the three forms 

fi{x, y, z) = {y, z, Ml + Ahz + Mgy + Bx + q{y, z)) , 

h{x, y, z) = {xo + -^x, yo - Miz, zq + M2Z + By + q{x, z)) , ^^^^^^ 

fsi^, y, z) = {xo + y + Mix, yo - -^x, zq + BM2Z + q{x, y)) , 

IVI2 

where the parameters have been rescaled. 

The dynamics of /2 and /s are trivial. Indeed for /2 the x dynamics decouples from the 
remaining dynamics, and for /s the (x, y) dynamics is linear and decouples from the z dynamics. 

Thus the only dynamically interesting case is fi. Since this map is the same as (jSJ, we have 
proved Th. EJ 

An example of this analysis consider the ACT-map (Arneodo-Coullet-Tresser map) which 
was introduced by Arneodo, Coullet and Tresser in their study of strange attractors in a family 
of differential equations on with homoclinic points of Shilnikov type |Du85j . 

X = ax — b{y — z), 

y = bx + a{y-z), (44) 
z = cx — dx^ + ez, 

where a,b,c,d,e are real parameters with J = e(a^ + 6^) / 0, 6d 7^ 0, and A; > 1 an integer. 
The ACT-map always has a fixed point at the origin. Nevertheless, when A: = 2 it is a quadratic 
family with a quadratic inverse, so it follows from Th.|l]that it is affinely conjugate to Q with 
Ml = 0. 

To see this explicitly define new coordinates using the transformation 

ax + by 

X = , Y = X , Z = ax - b[y - z) , 



a^+b"^ 



then, the map is rewritten as 



X = Y, 
Y = Z, 

Z = BX + SiY + S2Z - bdY^ 

where where B = e{a^ + 6^), = — (a^ + 6^ + c5 + 2ae), S2 = 2a + e. For the case k = 2 we can 
rescale the coordinates by a factor (bd)^^ to set the coefficient of Y'^ to one. Thus the map is 
written in the promised form @. 

Note that if we do one additional transformation and shift the variables by Si/ 2 we can also 
obtain the standard form @. Here we define Mi = jSi{2B + 2^2 + Si — 2), and M2 = S2- The 
inequality 

A = {l-B- M2f + 4Mi = {B + Si + S2-lf >Q 
holds guaranteeing the existence of a fixed point for the map. 
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8 Conclusions 



We have seen that many phenomena in the dynamics of three-dimensional Henon maps are 
markedly different from the possible phenomena predicted by our knowledge of two-dimensional 
dynamics. Primarily, we have seen these new phenomena when the Jacobian, B, is close to ±1. 
Two natural problems appear here, regarding the understanding of bifurcation scenarios that 
lead to chaotic three-dimensional dynamics. The first is to understand "how 2D dynamics is 
transformed into 3D dynamics" upon the variation of B from values near zero to those close to 
±1. The second problem — which is, in a sense, the inverse of the first one — is to understand "how 
3D volume-preserving dynamics is transformed into dissipative (and chaotic) dynamics" as B 
moves away from ±1. In the current paper (and also in jGOSTOSj), some partial results related 
to the second problem have been obtained; however, the first problem is — as of yet — unstudied. 

However, our analysis implies that these global problems have no unified, qualitative solution. 
One general conclusion is that the stable dynamics of the Henon maps of the first and second 
kind, and Q, have unmistakably different characters. The bifurcations in and ((21) both 
lead eventually to the formation of three-dimensional versions of Smale horseshoes. However, 
these have different characters: the horseshoes of contain type (2, 1) saddles of, while those 
of ^ contain type (1, 2) saddles. Thus, a third global problem appears: to understand "how the 
character of the attractors changes as one changes type of nonlinearity in 3D Henon families." 
In this case, it seems that using the family © with the parameters (a, b, c) determining the type 
of nonlinearity may be very profitable. 

While the problems posed above are of general interest, there are also many partial problems 
relating directly to the dynamics of quadratic maps. Some of these may also be solved in the 
the near future: 

1. Does the Henon map of the second kind, ((2j, have a large scale, wild- hyperbolic attractor? 

2. What are the dynamical properties of the local normal form for (j2|) near fixed points with 
three unit multipliers? 

This is related to the "reverse-Shimitsu-Marioka" system since ((2) is the inverse of 
Thus we will obtain the model (|33|) . with t replaced by —t. Consequently, the question is 
to determine the properties of the stable dynamics of this system. 

3. What are main bifurcations scenarios the 3D Henon map (pQ) that lead from simple dy- 
namics to a wild-hyperbolic, Lorenz attractor? 

There are two distinct, but useful approaches: (a) to fix B and consider the corresponding 
two-parameter family; (b) to consider B as the main governing parameter and to move for 
given values of (Mi, M2) from small values of B to values close to ±1. 

4. What role is played by the codimension-two bifurcations (recall ^ in the formation of 
attracting invariant sets? 

The analysis of codimension-one bifurcations alone gives us only a representation of the 
main hyperbolic elements of the dynamics. A detail analysis of codimension-three bifurca- 
tions (as in ^ allows us to find in local normal forms (using the fiow embedding approach) 
and verify the existence of small, wild-hyperbolic attractors. Analogous methods can be 
applied to the study of codimension-two bifurcations. We believe that this approach will 
provide an understanding of principal mechanisms for the creation of wild-hyperbolic at- 
tractors. 
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In summary, three-dimensional Henon maps have an abundance of new dynamical behavior 
which, it seems, has an important place in the theory of multidimensional dynamics as whole. 
We believe that they form a natural direction for a program of investigation similar to the well 
developed programs in one and two-dimensional dynamics. This is especially true since these 
maps cover (up to the affine conjugacy) all possible quadratic three-dimensional automorphisms. 
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